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A phase diagram is drawn in a parameter space of the nearly half-filled single band two- 
dimensional Hubbard model with U/t, t'/t and n as the parameters, U, t, t! and n being the 
on-site interaction, the nearest and second nearest neighbor transfer matrices and the number of 
electrons per site, respectively. The parameter space is chosen so as to include the models for the 
high-T c cuprates as well as 2D-organic superconductors. The superconducting and antiferromag- 
netic instability surfaces are drawn from the results of calculations by using the spin fluctuation 
theory within the fluctuation exchange approximation. The metal-insulator transition surface is 
also indicated in the diagram. From this phase diagram we discuss the possibility of the spin 
fluctuation mechanism being the common origin of the superconductivity in the high-T c cuprates 
and 2D-organic compounds. 

KEYWORDS: superconductivity, high-T c cuprates, organic compounds, spin fluctuations, Hubbard model, pseudo- 
spin gap 
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The origin of superconductivity in strongly correlated 
electron systems has been the central subject of current 
interest in condensed matter physics. It is widely be- 
lieved that the spin fluctuation mechanism is responsi- 
ble for the superconductivity in the heavy electron sys- 
tems. On the other hand, the mechanism for the high-T c 
cuprates is still under controversy. Although the mag- 
netic origin seems to be widely accepted, a long-standing 
controversy stilUaeisists between the spin fluctuation 
(SF) mechanisrroDff and the resonating valence bond 
(RVB) mechanism based on the t-J modeLo&fl'B 

Recently, the superconductivity in two-dimensional 
(2D) organic conductors: k-(ET) 2 A" [ET = BEDT-TTF, 
X = Cu{N(CN) 2 }A', X' = CI, Br, ...] has been success- 
fully. imzEslieated in terms of the spin fluctuation the- 
ory.0'BI___3E_l[!_P The conducting electrons in these sys- 
tems are believed to be well described in terms of a half- 
filled 2D-single band Hubbard model consisting of anti- 
bonding dimer orbitals. Thus the spin fluctuation mech- 
anism seems to be the only available mechanism which 
explainsJjbp anisotropic superconductivity in these sub- 
stances ;E3t3'E!f the t-J model as applied to a half-filled 
band leads to nothing but an insulating state. 

As for the cuprates, investigations in this decade have 
indicated that the essential physics is well described in 
terms of the 2D-single band Hubbard model with the 
transfer matrices up to the second or third nearest neigh- 
bors in the Cu02 plane. The more realistic d-p model 
consisting of Cu-d(x 2 — y 2 ) and O-pa orbitals is consid- 
ered to lead to essentially the same results as the Hub- 
bard model, or to be reduced, after proper renormal- 
izations of high frequency components, to an effective 
single band Hubbard modelEr It should be noted that 





Fig. 1. Transfer matrices in the 2D-Hubbard model (a) for 
cuprates and (b) for re-(ET)2-X", 2D-organic superconductors. 



the t-J model is derived from these models as the low- 
est order expansion from the strong coupling limit. Since 
the above-mentioned organic compounds are also well de- 
scribed in terms of the Hubbard model, we expect that 
these two classes of substances can be put in the same 
phase diagram by a proper choice of the parameter space 
for the 2D-Hubbard model. 

Let us consider a square lattice for both cuprates and 
organic compounds. The nearest neighbor transfer ma- 
trix is given by —t (t > 0). The second nearest neighbor- 
transfer matrix t' for cuprates as shown in Fig. 1(a) takes 
various values; the commonly used values are (t'/t) w 0.2 
for Lar_j.Sr2.CuO4 and ~ 0.5 for YBa 2 Cu 3 6+a ;. For 
brevity we neglect here the transfer matrices beyond the 
second neighbors. For k-(ET)2X the simplest and good 
model is to take only two out of the four second neighbor 
transfers t', as shown in Fig. 1(b). The estimated value 
for this transfer matrix is (t'/t) fss —0.8. The occupation 
number n of electrons per site is 1, or the band is half- 
filled, for k-(ET) 2 A^ and < 1 for hole-doped cuprates, 
best studied high-T c substances. So it seems illuminat- 
ing to draw a phase diagram in a three dimensional pa- 



1 



2 



Hisashi Kondo and Toru Moriya 



rameter space with n and t'/t for the horizontal axes and 
U/t for the vertical axis, taking the transfer matrices in 
Figs. 1(a) and 1(b) for t'/t > and t' /t < 0, respec- 
tively. For this purpose we must calculate the phase 
diagrams in various two-dimensional planes, U/t against 
n for various fixed values of t'/t, and/or U/t against t'/t 
for various fixed values of n. 

For this purpose, the information available on the 
ground state of the Hubbard model is extremely limited. 
Let us discuss the phase diagram calculated on the basis 
of the self-consistent theory of coupled modes of spin fluc- 
tuations as applied to magnetism and spin fluctuation- 
induced superconductivity. In practice, we use the sim- 
plest practicable one, i.e., the s^pJWl—Suctuatipn ex- 
change (FLEX) approximationHy-BBByM This 
approximation with no vertex correction has been suc- 
cessfully applied to the cuprates and the 2D-organic su- 
perconductors and recent studies on the vertex correc- 
tions resffijaL ihat they do not seriously affect the FLEX 
results.EJ'EaES' Since the numerical work of this type of 
theory is usually carried out at finite temperatures, we 
need to extrapolate the results to zero temperature. 

In our previous study of k-(ET)2X we showed such a 
phase diagram for U/t against the negative axis of t'/t 
for n = 1; the phase boundary between the paramagnetic 
metallic (PM) and superconducting (SC) phases was ob- 
tained by extrapolating T c against U/t curves to T c = 
for various ijalues of t'/t with the use of the Pade ap- 
proximants.EP A similar extrapolation of the T c against 
n curve to T c — for various values of U/t leads to the 
PM-SC phase boundary in the U/t against n plane for a 
given value of t'/t. The antiferromagnetic (AF) instabil- 
ity line or the PM-AFM (AF metal) phase boundary was 
obtained by extrapolating the inversed staggered suscep- 
tibility 1/xq to T = and picking up the parameter 
values giving 1/xq{T = 0) = 0. 

In order to construct a three-dimensional phase dia- 
gram as discussed above we must extend this type of 
calculation to include the models with various values of 
t'/t and n. The phase diagram for t'/t > is smoothly 
connected with that for t'/t < 0, with only half the sec- 
ond neighbor transfers, through the t'/t = plane. The 
former includes the cuprates and the latter includes k- 
(ET^A. We have calculated several phase diagrams as 
shown in Fig. 2, enabling the construction of a three- 
dimensional phase diagram, as sketched in Fig. 3. In 
Figs. 2 and 3 we adopt 10 > U/t > for the range 
of parameter values which is expected to cover both the 
cuprates and K-(ET) 2 A. 

Although this sketch is by no means accurate, the es- 
sential features of the phase diagram are considered to be 
reproduced correctly. The PM-SC phase boundary sur- 
face is a real one and should actually be observed. On 
the other hand, the antiferromagnetic instability surface, 
separating the PM and AFM phases is realistic when it 
is outside the PM-SC boundary surface, as is seen in the 
lower part of Fig. 3, but unrealistic when it is inside it, as 
is seen in the upper part of Fig. 3. In the latter case the 



3 



5 



3 



3 




7.5 
5.0 
2.5 

°-H 

10.0 
7.5 
5.0 
2.5 
0.0 



: 1 
• 




1 1 
* 


• 




• 


- '• SC 

AF 




sc y - 


_ PM 


ji....- 


PM _ 


(c)f'/f=0.0 \ 
i ■ 


i 


i 



0.8 



1.0 1.2 
n 



1.4 



i s 
♦ if 


^ J i 




• 


• i 






1 SC /' 


V. 




PM 


► 


(d) t'/t = 0.5 
i 


1 



0.6 o.; 



1.0 1.2 
n 



1.4 



Fig. 2. Superconducting and antiferromagnetic instability lines 
in (a) U/t against t'/t plane for n = 1 and in (b), (c) and (d) 
U/t against n plane for t'/t = —0.8, and 0.5. 



SC phase should extend beyond the PM-AFM boundary, 
since the SC phase is more stable than the AFM phase 
around this boundary. Thus we should find an SC-AF 
boundary surface, probably of first order transition. This 
surface, which is still to be calculated by comparing the 
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Fig. 3. Phase diagram in a three-dimensional (U /t)-{t' /t)-n 
space. Antiferromagnetic metal to insulator transition surface 
which should depend on the impurity potential is sketched here 
to give a general idea of its location. Phase boundary surfaces 
are cut at the plane of U/t = 10. 



energies of the two phases, is not shown in Fig. 3. The 
possibility of SC-AF coexistence is left for future inves- 
tigations. The antiferromagnetic insulator (AFI) phase 
exist in the n = 1 plane above an AFM- AFI boundary 
line located usually above the PM-AFM boundary line. 
In reality, the AFI phase extends into the 1 region 
owing to the Anderson localization due to the unavoid- 
able impurity potential. Considering this effect an AFM- 
AFI boundary surface is sketched arbitrarily in Fig. 3 as 
an aid to understanding the overall picture. 

From this figure we find that there is a substantial 
volume of superconducting phase outside the surface of 
antiferromagnetic instability in a range of intermediate 
size of U/t or U/W ~ 1, W « 8* being the bandwidth. 
As was mentioned previously, the superconducting phase 
for n = 1 (half-filled) can be understood only in terms of 
the spin fluctuation mechanism and this phase on both 
sides of the t'/t — plane smoothly continues to the 
electron-doped (n > 1) and hole-doped (n < 1) phases. 
The /c-(ET)2 compounds as well as high-T c cuprates find 
their positions within the one connected region of this 
phase diagram. This seems to indicate very strongly 
that these superconductors are consistently explained in 
terms of the spin fluctuation mechanism as was shown 
within the FLEX approximation. 

We have seen that the SC phase extends into the AFM 
region. It is interesting to note that the pseudo-gap phe- 
nomena, both in the spin excitation spectrum and in the 
one electron spectral density, are observed in this regime 



both in cuprates and in k-(ET) 2 XB»BBEZP Since these 
phenomena do not seem to be described by the FLEX ap- 
proximation we must improve the theory in this regime. 
Leaving discussions on the possible origins of the pseudo- 
gap phenomena for later part of this article, we may 
conclude for the moment that the SC phase outside this 
pseudo-gap regime is due to the antiferromagnetic spin 
fluctuations and that this SC phase smoothly extends 
into the pseudo-gap regime. 

Let us next consider to what extent the t-J model 
can cover this phase diagram. As mentioned previously, 
this model can never deal with the superconductivity for 
n = 1 as in k-(ET) 2 compounds. It is also hard to deal 
with the SC phase outside the AFM regime of relatively 
small U/t, at least below the critical value (U/t) c giving 
rise to an AFI phase. This part of the SC phase is well 
described by the spin fluctuation theory. What is clear 
is that the t-J model has a possibility of covering only 
a limited range of relatively large U ft and small 8n — 
|n — 1| in the phase diagram. 

The t-J model assumes an infinite value of U ft, strictly 
forbidding double occupancy of a site. Although a larger 
value of J corresponds to a smaller value of U/t, the 
complete exclusion of double occupancy makes it hard 
to judge if one can really approach a realistic intermedi- 
ate range of values for U / 1. As a matter of fact a typical 
value of J/t for the t-J model for the cuprates seems 
to be ~ 0.5, corresponding to U/t — 8 in the Hubbard 
model of the intermediate coupling regime. According 
to the variational Monte Carlo calculations for the su- 
perconducting condensation energy of these models, the 
Hubbard model with U/t ~ 8, t' = 0, n = 0.84 gives a 
value consistent with experiments, while the correspond- 
ing t-J model with J/t ~ 0.5 gives a value^ore than 
twenty times larger than the observed value. This re- 
sult seems to pose a serious problem for the t-J model 
which is an approximate form of the Hubbard model for 
large U/t and small Sn = \n — 1|. It seems that the part 
of the parameter space covered by the t-J model does 
not extend sufficiently far to overlap the part studied 
here which is expected to contain the high-T c cuprates 
as well as 2D-organic superconductors. 

Now let us discuss the pseudo-spin gap phenomena 
occurring in the transition region between the anti- 
ferromagnetic insulator phase and the superconducting 
phase. EIP Although there are several scenarios presented 
to explain these phenomena, none of them seems to 
be completely satisfactory. Recent observations of var- 
ious pseudo-gap behaviors in K.-tVT^pX quite similar 
to those in underdoped cupratesaEJ'E^LljEP indicate 
strongly that the origin of the phenomena is common 
to both of these systems. 

We stress here that the possible mechanism to explain 
the pseudo-gaps in k-(ET)2^ is naturally expected to 
cover those in cuprates, too. On the other hand, the 
t-J model mechanism or the spinon pairing mechanism 
for the pseudo-gaps cannot cover k-(ET)2X . Thus it 
seems wise not to pursue the t-J model scenario for the 
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pseudo-gaps before the range of validity of the model 
itself is clarified. In doing so we can avoid looking for 
two different mechanisms for the same phenomena. 

The other proposed mechanisms include the spin den- 
sity wave (SDW) fluctuations, the coexistence of charge 
and spin density wave (CDW and SDW) fluctuations and 
the influence of superconducting fluctuations or forma- 
tion of preformed pairs above T c . We note that these 
mechanisms for the pseudo-gaps, the validly of which is 
still to be examined, necessarily assume the spin fluctu- 
ation mechanism as the origin of the d-wave supercon- 
ductivity. 

As a consequence of the 2D-antiferromagnetic spin 
fluctuation theory we expect to have 2D-SDW fluctu- 
ations of very low frequencies beyond the AF instability 
point. This seems to give a natural explanation for the 
pseudo-gap in one-electron spectral density. A difficulty 
with this mechanism is that it is incapable of explaining 
magnetic excitation gaps or the temperature dependence 
of the nuclear spin-lattice relaxation rate. Insofar as we 
have an AFM or SDW ground state (when the SC phase 
is suppressed) this difficulty seems unavoidable. One so- 
lution may be to consider additional charge density wave 
(CDW) fluctuations which provide a magnetic excitation 
gap. Although this is an appealing possibility for both 
cuprates and organic compounds, it is not clear if the 
CDW fluctuations can be so universal as the pseudo-gap 
phenomena in these systems. 

We may think of another possible way out of this dif- 
ficulty. Unfortunately the AFM phase in a single band 
2D-Hubbard model has not been well enough studied 
beyond the mean field and FLEX approximations. Even 
the stability of the metallic AF ground state in this case 
has not been proved either theoretically or experimen- 
tally. Possibly the necessarily small AF moment (< 1/2) 
is not sustained in this 2D-metallic phase with large 
quantum fluctuations, leading to a non-magnetic ground 
state with a pseudo SDW and magnetic excitation gaps. 
The pseudo-SDW gap formation would naturally sup- 
press the superconducting gap, leading to an explanation 
for the doping concentration dependence of T c in the un- 
derdoped regime. The key issue here is to explain the 
magnetic excitation gap, possibly by local singlet corre- 
lations or magnetic anisotropy. Although this is simply 
conjecture to be substantiated later, the influence of this 
possible mechanism on the spin fluctuation-induced su- 
perconductivity will not be large since high frequency 
components of the spin .fluctuations mainly contribute 
to the superconductivityEil'Lj' 

Lastly, as a consequence of the above arguments we 
suggest that it seems unlikely that the t-J model expla- 
nation for the high-T c superconductivity in cuprates and 
the spinon pairing mechanism for the spin pseudo-gap 
phenomena are realistic. Further careful examination of 
the range of validity of the t-J model is desired. 

We may now conclude with high probability that the 
superconductivity of the organic 2D-compounds as well 
as that of the high-T c cuprates is consistently explained 



in terms of the spin fluctuation mechanism. Although we 
still have one serious problem of explaining the pseudo- 
gap phenomena observed in the underdoped cuprates 
and in the organic 2D-compounds, the main physics of 
the spin fluctuation-induced superconductivity is not ex- 
pected to be much influenced by the underlying mecha- 
nism of these phenomena. 
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